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This paper describes a complete analytical formalism for calculating electron subband energy 
and degeneracy in strained multi-valley quantum wells grown along any orientation with explicit 
results for AlAs quantum wells. In analogy to the spin index, the valley degree of freedom is 
justified as a pseudospin index due to the vanishing intervalley exchange integral. A standardized 
coordinate transformation matrix is defined to transform between the conventional-cubic-cell basis 
and the quantum well transport basis whereby effective mass tensors, valley vectors, strain matrices, 
anisotropic strain ratios, piezoelectric fields, and scattering vectors are all defined in their respective 
bases. The specific cases of (001)-, (110)-, and (lll)-oriented aluminum arsenide (AlAs) quantum 
wells are examined, as is the unconventional (411) facet, which is of particular importance in AlAs 
literature. Calculations of electron confinement and strain in the (001), (110), and (411) facets 
determine the critical well width for crossover from double- to single-valley degeneracy in each 
system. The biaxial Poisson ratio is calculated for the high-symmetry lower Miller index (001)-, 
(110)-, and (lll)-oriented QWs. An additional shear strain component arises in the higher Miller 
index (411)-oriented QWs and we define and solve for a shear-to-biaxial strain ratio. The notation 
is generalized to address non-Miller-indexed planes so that miscut substrates can also be treated, 
and the treatment can be adapted to other multi-valley biaxially strained systems. To help classify 
anisotropic intervalley scattering, a valley scattering primitive unit cell is defined in momentum 
space which allows one to distinguish purely in-plane momentum scattering events from those that 
require an out-of-plane momentum component. 

PACS numbers: 73.21. Fg,73.50.Bk,73.61.Ey,71.70.Fk,71.70.Gm 



I. INTRODUCTION 



Calculating valley degeneracy in a quantum well re- 
quires a comprehensive treatment of strain, quantum 
confinement, and piezoelectric fields since all contribute 
at comparable energy scales. Of the various multi-valley 
semiconductors, the indirect-bandgap zinc blende semi- 
conductor aluminum arsenide (AlAs) with its bulk three- 
fold valley degeneracy is of particular interest (Fig. [U be- 
cause its heavy anisotropic electron mass allows for large 
interaction effects^, and its near perfect lattice match 
to GaAs substrates allows for high-mobility, modulation- 
doped quantum wells (QWs) 2 .^. AlAs/AlGaAs QWs can 
reach mobilities of the order of // = 100,000 cm 2 /Vs 
in (OOl)-facet QWs^ and in the high mobility direc- 
tion of anisotropic (llO)-facet QWs^. Unconventional 
facets such as (411) have also proven useful in identify- 
ing exchange effects like quantum Hall ferromagnetism 
in AlAs QWs 6,7 . Evidence for a QW width crossover 
from double- to single- valley occupation has been shown 
for (OOl)AlAs wells 8 ^, as has evidence for single valley 
occupancy in wide (110) AlAs wells^. Dynamic control 
of the valley degeneracy has been realized with uniaxi- 
ally strained (OOl)AlAs QWs to induce valley degener- 
acy splitting ^ 10 ' 11 . Such studies can quantify interac- 
tion effects, calibrating valley strain susceptibility and 
valley effective mas a 10 i 12 ~ — . Quantum confinement to a 
one-dimensional multi-valley system has been achieved 
in cleaved-edge overgrown quantum wire o 15 ' 16 and quan- 
tum point contacts^ 7 .. Novel interaction effects in QWs 
include anisotropic composite Fermion mass 18 and valley 



skyrmions, whereby electrons populate linear superposi- 
tions of two valleys at onceiS. Such interaction effects 
that result from exchange splitting of a perfect SU(2) 
symmetry 2 ^ may prove useful in future quantum device 
applications, where the valley degree of freedom functions 
as a pseudospin. 

Various formalisms have been developed to understand 
valley occupancy in a multi- valley system. AlAs QWs are 
slightly strained with respect to the GaAs substrate, so 
one must consider both quantum confinement and strain 
to estimate the electron subband energy and degeneracy 
of the different electron valleys, and certain growth sym- 





FIG. 1: (Left) Real space depiction of AlAs QW structures 
defining the transport basis axes, a, 6, c. (Right) Brillouin 
zone for bulk AlAs defining the conventional-cubic-cell ba- 
sis x,y,z. Note the three degenerate valleys that are oc- 
cupied with electrons as indicated by ellipsoidal equi-energy 
contours. 
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metries also introduce in-plane or out-of-plane piezoelec- 
tric fields. Stern and Howard modeled two-dimensional 
(2D) confinement of electron valleys in Si for arbitrary 
crystal plane orientations neglecting strain effects 21 . Van 
de Walle theoretically studied the absolute energy level 
for an unstrained semiconductor heterojunction, con- 
sidering strain effects only for bulk systems 2 ^. Smith 
et al. calculated the strain tensor theoretically for the 
case of (001)- and (lll)-oriented QW superlattices^. Al- 
though Caridi and Stark derived the complete strain ten- 
sor for arbitrarily oriented substrates with cubic symme- 
try, they ignored a critical shear component which arises 
in the high Miller index directions under relaxation 24 . 
De Caro et al. calculated the shear strain with com- 
mensurability constraint equations accurately only for 
the low index (100)-, (110)- and (lll)-orientations 25 . 
Yang et al. corrected the constraint equations and cal- 
culated piezoelectric fields arising from the shear strain 
components for the general case of a pseudomorphic 
film 2 ^. Adachi calculated the strain tensor theoretically 
for the case of bulk and superlattice structures but did 
not address the case of single QW structures 2 ^. Ham- 
merschmidt et al. focus on the case of isotropic planar 
strain arising in QWs but do not carry out a combined 
treatment of strain energy along with quantum confine- 
ment energy and how these affect valley degeneracies 2 ^. 
Theoretically, Rasolt proposed an understanding of the 
various continuous symmetries and symmetry-breaking 
in a multi-valley system in terms of an SU(iV) symme- 
try where iV is the valley degeneracy number—, but did 
not explicitly derive the vanishing intervalley exchange 
integral which underlies this theory. Material specific 
calculations of quantum confinement and strain in high- 
mobility Si and Ge structures are well studied for high- 
symmetry facets 2 ^— and recently, such calculations have 
been been made in other quantum confined systems as 
wel l 41 ' 42 . However, a combined treatment of quantum 
confinement, strain, and piezoelectric fields to determine 
valley degeneracy in QWs is lacking, especially for low 
symmetry facets. 

To eventually model transport in such valley degener- 
ate systems, one must also consider the extra scattering 
channel not present in single valley systems, namely in- 
tervalley scattering. It is therefore useful to introduce a 
/c-space unit cell which permits visualization of momen- 
tum scattering events in a geometry that is natural to 
the quantum confinement direction. However, the stan- 
dard depiction of a 2D Brillouin zone of a quantum con- 
fined valley-degenerate system^!, projects all valleys to a 
single two-dimensional plane. Such a depiction loses in- 
formation about the out-of-plane momentum scattering 
component that was projected out, which is necessary to 
determine the full momentum scattering matrix element. 
Thus it is useful to develop a graphical representation 
for the unit cell in fc-space that can clearly elucidate how 
valleys are coupled with both in-plane and out-of-plane 
momentum-scattering events. 

The paper is organized as follows. In Section [Hi ex- 



change energy calculations demonstrate why the valley 
index can function as a pseudospin index. In Section UTT1 
we develop the notation and formalism for determining 
valley degeneracy in multi- valley strained semiconductor 
QWs. The three subband energy components are defined 
- kinetic energy, confinement energy and strain energy - 
as are the two useful coordinate bases - the conventional- 
cubic-cell basis and the transport basis - as well as the 
coordinate transformation matrix which transforms be- 
tween them. Single-electron analytical solutions are pro- 
vided to allow easy identification of characteristic en- 
ergy scales for multi-valley systems. The explicit crite- 
ria for valley degeneracy in arbitrarily oriented biaxially 
strained QWs are defined in Section [TV] Section [V] de- 
scribes how shear strain can induce piezoelectric fields in 
the QW. In Section lVTl we introduce the valley-scattering 
unit cell - a 3D primitive unit cell with the same volume 
as the Brillouin zone, but which shares the full symme- 
try of the reciprocal lattice vectors that lie in the Miller 
index plane. Momentum-space illustrations provide in- 
tuition for visualizing intervalley scattering. Section IVIII 
applies the developed formalism to the case of AlAs QWs, 
whereby the projected in-plane transport masses, strain 
tensor, degeneracies, piezoelectric fields, and valley split- 
ting energies for the specific cases of (001)-, (110)-, (111)- 
as well as the unconventional (411)-oriented QWs illus- 
trate how valley degeneracies can be engineered. We re- 
view the crossover width calculation for double-to-single 
valley occupation in the (001)-, (110)-, and (411)-facet 
QWs. For the high-symmetry low Miller index orien- 
tations (001), (110) and (111) we calculate the biaxial 
Poisson ratio, and for the low-symmetry high Miller- 
index (411) facet there arises a shear component in the 
strain tensor and consequently a shear-to-biaxial strain 
ratio. We calculate intervalley scattering ratios for valley- 
degenerate AlAs QWs on various facets in Section [Villi 
and refer to the valley scattering unit cell for intuition. 
We complete our analysis by comparing the standard 2D 
Bravais lattice valley representation to our valley cell rep- 
resentation. To generalize the formalism for calculating 
valley subband energies in arbitrarily oriented facets, we 
end in Section 11X1 by adapting this notation for miscut 
samples 4 ^. 



II. VALLEY EXCHANGE INTEGRAL AND 
VALLEY INDEX AS A PEUDOSPIN 

Just as electrons with different spins are distinguish- 
able and have no exchange term in their interaction en- 
ergy, electrons in different valleys can be shown to be 
effectively distinguishable from each other due to the van- 
ishingly small intervalley exchange energy. The valley in- 
dex can thus be treated as a pesudospin index. Rasolt 20 
discusses symmetry-breaking in a multivalley system but 
does not explicitly derive this exchange integral. Previ- 
ous work on Si quantum dots to study Kondo effects- 
calculates the exchange integral only for the special case 
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of identical spatial wavefunctions on Si quantum dots. 
For completeness, in this section we explicitly derive the 
inter- and intravalley exchange integrals for arbitrary 
wavefunctions to justify the valley index as a pseudospin, 
with exchange interaction proportional to 8 t y where r 
and t are valley indices. 

An electron wavepacket within a single valley is de- 
scribed with a weighted integral over Bloch functions. 
The normalized wavepacket -0 TjCr in valley r is chosen to 
be in a spin polarized state a described by the spinor \ a 
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J d 3 k A k e ^ T + k ) r u q . +k (r) Xtr . (1) 



M q r +k (r) is the component of the Bloch function peri- 
odic in the Bravais lattice, A-^ is the complex amplitude 
for a particular k, q T + k is the total crystal momentum 
where q T denotes the center of the valley and k denotes 
the additional small momentum deviation away from this 
valley center, and V is the volume of the system. Assum- 
ing k is much smaller than the Umklapp vector for the 
lattice, and taking the envelope function approximation 
Wq T +k — u q r, we obtain 



UqT (r)e^- r X ^(r) 



(2) 



where <fi(r) is the slowly varying envelope function. The 
Coulomb exchange energy integral between the two elec- 
tron wavefunctions ^ T)(7 (r) and 2p- w (r) in valleys r and 
r with coordinates ri and r 2 is given by 



= E / / dr ^ dr 



[ -2] rVv,<x(ri)VC,<T( r 2) 

ri - r 2 



^r,a( r 2)Vv,a(ri) 



(3) 



Substituting Eq. ^ in Eq. ©, we can show that two 
particles with different valley indices behave like distin- 
guishable particles by virtue of their vanishing exchange 
integral, 



E. 



T.T.O.O U (T(T 



V 2 



l r i - r 2| 



= i Q TT -(r 2 -r 



dv x dr 2 cj){r 1 )(j)*{v 2 )(j){Y2)cj) (ri) 
l} u qT (ri)w* (r 2 )u q -(r 2 )u*-(ri) 



with intervalley scattering wavevector 
Q T " = q r q T , 



(4) 



(5) 



where the inner product of the spinors is given by XwX<? ~ 
5aa- This delta function denotes that the exchange in- 
tegral vanishes when the two spin indices are different 
u/F. Analogously, if r ^ f, Q TT is of order an Umk- 
lapp vector, and the rapidly oscillating complex expo- 
nential makes the integral vanish for envelope functions 



larger than a few lattice constants. On the other hand if 
t = r, Q TT = and the exponential term is unity and 
the integral remains finite. Therefore, the dependence 
of the exchange integral Eq. (Q} on valley index can be 
approximated with a second delta function to notate the 
vanishing exchange integral between different valleys 



V" 2 



dr 2 0*(r 2 )0(r 2 )- 



V r 2 



(6) 



We note that Eq. ^ cannot be simplified because the 
Coulomb potential can vary rapidly over small distances 
of order a lattice constant. 

We conclude by virtue of this vanishing intervalley ex- 
change interaction, that wavefunctions in different valleys 
can be treated like distinguishable particles, and the val- 
ley index is a valid pseudospin index. As was pointed out 
by Rasolt 20 the valley pseudospin constitutes an SU(iV) 
group where N is the number of electron valleys. 



III. VALLEY SUBBAND ENERGY 

Crystal symmetry dictates that multiple energy- 
degenerate valleys will occur whenever a local conduc- 
tion band minimum exists away from the origin in mo- 
mentum space (the T-point). When such a multi- valley 
system is then quantum confined in a layer with planar 
Miller indices M = (hkl), the energy E T (k) in the lowest 
subband in valley r is 



E T (k) = £ r (k)+T r (k) + A£ T 



(7) 



where k is the 2D in-plane momentum relative to the 
t- valley minimum, -Eg(k) is the ground confinement en- 
ergy, T T (k) is the in-plane kinetic energy, and AE T is the 
strain induced energy shift caused by lattice mismatch of 
the QW with respect to the substrate. Note that the 
Miller index M is not explicitly superscripted because 
it is common to all valleys. To calculate these terms in 
Eq. 0, we need to find the in-plane (parallel to the QW) 
and out-of-plane (confinement direction) components of 
the inverse mass tensor of the corresponding electron val- 
ley as well as the various strain tensor components in the 
QW. 

We start by introducing two useful bases, the 
conventional-cubic-cell (CCC) basis x = (x, y, z) and the 
transport basis a = (a, 6, c), along with the coordinate 
transformation matrix R M which transforms between 
them 



a = R x 



(8) 



The CCC basis has the x, y, and z axes aligned along the 
axes of the cubic cell of the reciprocal lattice of the crys- 
tal, Fig.[T] (right). The transport basis is chosen with the 
a-b plane parallel to the QW, and the c-direction perpen- 
dicular to this plane, Fig. Q] (left). Thus if the Miller index 
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of the growth plane is (hkl), the perpendicular unit vec- 
tor in the transport basis is c = (h, k, I) v / fe2+fc2+f5 ■ To 
uniquely define the a and b directions in the transport 
basis, we take a to be the in-plane unit vector with the 
lowest Miller index, and b the unit vector which main- 
tains a right handed coordinate system a x b = c. With 
this definition we obtain a unique set of axes. These also 
define the components of the coordinate transformation 
matrix R M whereby a, b, and c are the top, middle, 
and bottom rows of the coordinate transformation ma- 
trix R M . In what follows, vectors and tensors which are 
expressed in the x-basis are unprimed, and in the a-basis 
they are primed. 

The mass tensor is naturally expressed in the unprimed 
CCC x-basis, so by transforming it to the primed trans- 
port a-basis we can easily extract the in-plane and out- 
of-plane confinement masses for a given electron valley. 
The matrix inverse of the mass tensor for the r-valley 
in the x-basis is denoted by w T = (m T ) _1 . We as- 
sume parabolic bands so that the mass is independent 
of the wavevector k. Applying the transformation (w' T ) 
= R M (w r )(R M ) _1 we obtain the inverse mass tensor of 
the r-valley in the primed transport a-basis 



(9) 
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2 sub-matrix w' £ X 2 gives the inverse 



The upper left 2 x 
in-plane mass tensor and w' T cc gives the inverse of the 
confinement mass perpendicular to the QW. 



A. Quantum confinement energy 

We can now define the first energetic term in Eq. (|7|), 
the quantum confinement energy E^k), using the co- 
ordinate transformation matrix and out-of-plane compo- 
nent of the inverse mass tensor in the primed a-basis. 
This is calculated from the ground state energy solution 
of the ID Schrodinger equation for a particle confined 
along c-direction with in-plane momentum k, 



-h 2 d 

2 dc 



w 



#k(g) 

dc 



(10) 

where w' cc (c) is the diagonal component of the recipro- 
cal mass tensor in the a-basis. Due to the different re- 
ciprocal masses in the barrier (w' g) and the well (w' 1^) 
layers, wavefunction derivatives at the boundary Co must 
satisfy^ 



dc 



IJ ' W,cc" 



dipk(c) 



dc 



(11) 



The confinement potential V(c, k) can be expressed as 



V(c,-k) = 



if |c| < W/2; 
V (k) if |c| > W/2 



(12) 



where the height of the potential barrier Vq is given by 

2,w)'k (13) 



V (k) = £ c T , B -£c T ,w+yk.(w'L 2 ,B-w'L 

where -E T W is the conduction band energy of the r-valley 
in the QW, E^ B is the energy in the barrier, and W is 
the well width. The last term accounts for differences 
in the in-plane effective mass as a k-dependent barrier 
height. This last term is often neglected in calculations 
under the assumption that k is small and the tensor mass 
difference w' B — w' ^ is also small, making Eq and ip( c ) 
independent of in-plane momentum. 



B. Kinetic energy 

We can now calculate the second term in Eq. ([7]), the 
kinetic energy T(k), from the in-plane inverse mass ten- 
The in-plane inverse mass tensor is deter- 



sor w 



2x2- 



mined by solving the ground state from Eq. (|10p and then 
determining a weighted average of the well and barrier 
inverse mass tensors given by 



W/2 

2X2 = 2w '2x2,W / \lj}(c)\ 2 dc + 



2w '2x2.B / dc 

W/2 



(14) 



In most cases it is safe to assume the wide- well limit, 
whereby the effective mass is equal to that in the well 
material only. The kinetic energy is 



T-(k) = T k.(w^ x2 ).k 



(15) 



The projected in-plane effective masses can be solved 
from the determinant of this same sub-matrix 



det(w'^ x2 - AI) = 



(16) 



The eigenvalue solutions Ai and A2 directly give the re- 
ciprocal masses along the major and minor axes of the 
projected mass tensor. The cyclotron mass is the same 
as the 2D density-of-states mass, and is given by 



'2D 



= det(w'J x2 )- 1 / 2 = (A^)- 1 /? 



(17) 



The spin-degenerate two-dimensional energy density-of- 



states per unit area in the r-valley is 



'2D 



that the total energy density-of-states is given by a sum 
over all valleys 



where Q(x) is the Heavyside step function. 



(18) 
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C. Strain energy 

We can now define the final term of Eq. ([7]) , the strain 
energy AE T , in terms of the strain tensor and the co- 
ordinate transformation matrix. Following the notation 
of Van de Wallet and Herring and Vogi>^, the absolute 
energy shift AE T of the r-valley for a homogeneous de- 
formation described by the strain tensor e in the x-basis 
is given by 



22.50 



AE T = (Z T d 6 i: 



y 



(19) 



where q r is a unit vector in the direction of the r-valley, 
and and SJJ represent the deformation potentials due 
to a bulk dilation and a uniaxial deformation, respec- 
tively. The average shift in the energy of the subband 
extrema is given by 



AE av = (SJ + l -~l)kz lJ = «c 



(20) 



where a c is the hydrostatic deformation potential for the 
conduction band and £ye y = Tr(e) is the trace of the 



strain tensor t xx + e.„ 



Often the relative strain 



energy shift A T from the mean value is all that is needed 
to determine subband occupancy: 



A T = AE T - AE a 



l 6 ^ 



(21) 



To determine the strain tensor e it is easiest to explic- 
itly determine its components in the a-basis, which we 
will denote with e', and then apply a rotational transfor- 
mation to express it in the x-basis for use in Eqs. f|19[)- 
([21]) . Heteroepitaxial QWs are biaxially strained relative 
to an unstrained substrate with a different lattice con- 
stant. Thus the strain tensor e' is diagonal in the upper 
2x2 block e' aa — e'bb = e'n, where e'ii is the in-plane 
strain component given b y 22 i 



^substrate ^laycr 
^laycr 



(22) 



Note that lattice matching to the unstrained substrate 
fixes e' ab = e' ha = 0. The remaining components of 



-6a 



the strain tensor are linearly proportional to the parallel 
strain e'„ as follows: 



M 



u 



(23) 



Using Hooke's law in the x-basis, stress 



constants 

(er) and strain (e) tensors are related by 

&ij = Cijki e k t i,j,k,le {x, y, z] 



(24) 



Here c is the fourth-rank elastic stiffness tensor in the 
x-basis. For crystals with cubic symmetry c can be sim- 
plified using the Voigt notation^ into a matrix form 



Oi = cv 



i,j = 1,2,..., 6 



(25) 



The indices 1 through 6 denote xx, yy, zz, yz/ zy, zx/xz 
and xy/yx, respectively. For cubic materials most of 
the elements of the matrix C r ; vanish and Cn = C22 = 



C337 C12 — C13 — C23 and C44 — C55 
the stress tensor er to 
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(26) 



where ej = 

and €6 = 2e xy . Table Q] lists the values of Cii,Ci2, and 
C44 for AlAs and GaAs. The strain tensor e and the 
elastic constants cumn are used to obtain the free energy 
of isothermal clastic deformations of a medium 2 ^, 



F{e) 



(27) 



klmn 



For a cubic crystal using the Voigt notation for the elastic 
constants, the elastic energy reduces to 



F(e) = ^f(e 2 xx + 4 + 4) + 2C 44 (4 + 4 + 4*) 

+ Cl2( ^xx^yy H~ ^xx^zz H - ^-yy^-zz ) (28) 



The coefficients Df 1 can be determined for a given 
crystal facet orientation by transforming the strain ten- 
sor of e' into the x-basis with the rotation matrix 



e = (R M )- 1 (e')R M 



(29) 



and inserting these components into the elastic energy of 
Eq. (|2"5)l and minimizing with respect to each Df 1 



F(e) =F[(TL M )- 1 (e')R M ] 



(30) 



where Dq 1 = — e 



/e' aa = — £y/£j_ is the biaxial Poisson 
ratio (notated D M by Van de Wallet), and we define 



M 



le' 



and Df = e'Je'a, 



as the shear-to-biaxial 



strain ratios. 

The coefficients Df 4 can be derived for a crystal with 
cubic symmetry and arbitrary Miller index M by mini- 
mizing the free energy of the layer in terms of its elastic 



dF(e) 



(31) 



to deduce a set of three equations which can be simul- 
taneously solved to give the Df values. With the Df 
values solved for this material and this facet, the strain 
tensor in the x-basis can now be used in Eqs. (|19l) - (|2"Tj) 
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material 


D° 01 




Co" 1 


Df 1 


Df 1 


[ 


Cl2 

GPa] 


C44 


a c 
[eV] 


[eV] 


mi 
[m e ] 


nit 
[m e ] 


e x ,4 
[C/m 2 ] 


AlAs 


0.854 


0.616 


0.550 


0.775 


0.176 


125.0 


53.4 


54.2 


2.54 


6.11 


1.1 


0.20 


-0.23 


GaAs 


0.934 


0.580 


0.489 


0.820 


0.250 


122.1 


56.6 


60.0 


-0.16 


8.61 


1.3 


0.23 


-0.16 



TABLE I: Strain ratios Df 1 for different facet orientations, elastic constants Cy, deformation potentials a c ,S^, X- valley mass 
tensor components mi, rat, and piezoelectric coefficients e x ,4 for AlAs and GaA a 21 ' 22 ' 52 ^ . 



to deduce the strain energy contribution to each valley 
AE T . 

As discussed in Section IVII1 for high-symmetry crys- 
tal orientations like (001), (110), and (111), there is no 
shear component in the strain tensor when expressed in 
the transport a-basis so that Df* = D^ 1 = 0, whereas 
for (411) growth the Dif coefficient is nonzero, and the 
direction of shear is defined by the vector a given by 

a = Df a + Dfb . (32) 

For AlAs and GaAs, we tabulate the nonzero Df 1 in 
Table U along with the elastic stiffness components in 
Voigt notation. For completeness, the Appendix provides 
a compact list of equations for determining the strain ten- 
sor for arbitrary substrate orientation. From the above 
analysis, we can now calculate all the terms in Eq. ([7]) 
for any valley and any QW orientation. 

IV. ROBUST VALLEY DEGENERACY AND 
CROSSOVER VALLEY DEGENERACY 

Two valleys r and f will be degenerate if 

£ T (k = 0) =E T (k = 0) . (33) 

When this condition holds for r and r regardless of the 
QW width, we call the valley degeneracy robust, and if it 
holds only for a specific well width Wo , we call it crossover 
degeneracy. The crossover degeneracy condition can arise 
only when the valley electrons with a higher confinement 
mass are strain-shifted to a higher energy than the valley 
electrons with a lower confinement mass. 

Robust valley degeneracy can only occur if the strain 
energies and out-of-plane confinement masses are inde- 
pendently equal. The strain energy must satisfy 

AE T = AE T , (34) 

for all well widths, which simplifies via Eq. (fTi?]) to the 
x-basis condition 

q- e q- = q-. e .q- (35) 

or equivalently in the a-basis, 

q' T • e' • q' T = q' T • e' • q' T . (36) 

When D^ 1 = D^ 1 = 0, this condition for degeneracy 
simply reduces to 

|q' " - c| = - c| , (37) 



meaning that all valleys with the same polar angle from 
the c-axis are degenerate. For the more general case of 
arbitrary Df 1 , the full Eq. (|3"o| has to be solved to deter- 
mine valley degeneracy. However, given the shear vector 
a in the plane of the QW as defined in Eq. (f3"2"j) . a spe- 
cial case can be defined, and two valleys r and r are 
degenerate if they simultaneously satisfy Eq. (f3"Tf and 

\q' T -a\ = |q'" ■ S| , (38) 

or equivalently, if the valleys are mirror symmetric about 
the c-a plane. 

The second condition for robust valley degeneracy is 
that the out-of-plane inverse masses are equal, 

w'lc = W 'L ■ (39) 

guaranteeing equal confinement energies. As long as the 
mass ellipsoid associated with each valley is oriented with 
its longitudinal mass parallel to the unit vector q r , as 
is usually the case, then the Eq. (f39|) is automatically 
satisfied under the same condition as Eq. (|37|) . 



V. PIEZOELECTRIC EFFECTS ON THE 
QUANTUM WELL POTENTIAL 

In crystals which are inversion asymmetric like zinc 
blende crystals, strain may also produce piezoelectric 
fields, modifying the confinement potential V(c, k), and 
resulting in shifts of the ground energies Eq for the var- 
ious valleys as solutions to the Schroedinger equation, 
Eq. (fTT)]) . The out-of-plane piezoelectric field which re- 
sults will not break robust valley degeneracy, but will 
shift the crossover degeneracy condition to a different 
well width, Wo. 

In our treatment, we will calculate the piezoelectric 
field inside the QW and assume that any out-of-plane 
component will be canceled with an external gate bias, 
restoring the condition of the flat square well. The piezo- 
electric field inside the QW and/or the barriers is propor- 
tional to the shear strain in the crystal basis^i 



Si = aj €j . (40) 

Here £ , lists the unprimed piezoelectric field components 
in the x-basis indexed by i — x, y, z, Eq and s s are 
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the free-space and relative semiconductor dielectric con- 
stants, and eij represents the piezoelectric tensor with j 
indices in Voigt notation. We assume that £ , is zero in the 
unstrained substrate. Recall that e = (R M )" 1 (e')R Jl/ 
is the strain tensor derived from Eq. ([23]) but expressed 
in the unprimed x-basis. Even though the facet orienta- 
tions (110) and (111) have no shear strain in the growth 
a-basis, they do have shear strain in the crystal x-basis. 
The (001) facet orientation, on the other hand does not 
have any shear strain. The strain in Eq. (|40p is repre- 
sented in Voigt notation as for Eq. ([21))) with j — 1, 2, ...6. 
For zinc blende crystals, the only nonzero piezoelectric 
coefficients are e x> 4 — e Vt s = e z ,(r^ and we obtain 



f 






fe 1 


\ 












£3 








£4 




xA / 






£5 






u 


J 



(41) 

Table Q] lists the piezoelectric coefficients for GaAs and 
AlAs, and the coefficient for Al^Gai-^As can be linearly 
interpolated according to the x-Al contend. The general 
shear strain components in a QW are reproduced from 
Cadini and Star h 24 i 55 and simultaneously derived using 
Eq. ([29]), Eq. ([23]), Eq. ([30]) and Eq. ([31]). It is useful to 
express the electric field in the primed growth a-basis 




-2e 7 



R 



A I 



t-xy 



(42) 



The c-component of the piezoelectric field £' c will affect 
the electron subband energy by modifying the confine- 
ment potential V(c, k) of Eq. (|12[) both in the strained 
well and in the strained barrier 



V(c,k) 



(43) 



Vb(k) 
I Vo(k) 



Q £ ' c ,w c 
, . w W/2 + q£' c>B (c- 



q£' c 



if |c| < W/2; 
W/2) if c> W/2; 
W/2) if c < -W/2, 



where £' c w and £' c B are the out-of-plane c-components 
of the electric field in the well and barrier, respectively. 
As pointed out by Adachi^l for (111) GaAs and AlAs, 
this piezoelectric field points from the (lll)B surface (As- 
terminated) to the (lll)A surface (Ga-terminated) when 
under positive c-axis strain. 

The piezoelectric field in the QW and barrier mate- 
rial alters both the quantum confinement energy and the 
wavefunction, and in the limit of weak piezoelectric fields 
the energy shift can be estimated to be Ep Z = £' c W/2. For 
example, in strained (111) AlAs on a GaAs substrate, the 
piezoelectric field is £' c = 8.794 x 10 6 V/m, so for a QW 
width of W = 10 nm the piezoelectric energy shift would 
be around = 43.97 meV. 



For the present treatment, we reiterate our assumption 
that any nonzero out-of-plane component of piezoelectric 
field £' c will be assumed to be canceled by an external gate 
voltage. 



VI. VALLEY SCATTERING UNIT CELL 

In addition to calculating valley subband energies and 
degeneracies, it is useful to visualize the valley orienta- 
tions in three dimensions in order to map relevant inter- 
valley scattering vectors. In this section, we lay down a 
general derivation of a 3D primitive unit cell which we 
call the valley scattering unit cell (VSC), that aids in 
identifying in-plane and out-of-plane intervalley scatter- 
ing vectors. The proposed VSC is not a 3D Brillouin 
zone, but a simple unit cell defined to have the 2D sym- 
metry of the Miller index plane of the QW, while pre- 
serving the total volume of the original 3D Brillouin zone. 
This cell is to be also distinguished from the standard de- 
piction of a 2D Brillouin zone which neglects information 
about the out-of-plane momentum which is necessary for 
calculating momentum scattering matrix elements. 

To define such a unit cell, we first identify the 2D sub- 
set of reciprocal Bravais lattice points which lie in the 
Miller index plane of the QW. The original 3D Bravais 
lattice is composed of parallel planes of this 2D sublat- 
tice which are displaced and separated by an interplanar 
wavevector 2ir/d. We now take the primitive unit cell of 
this 2D sublattice and extend it in the perpendicular di- 
rection by an amount ±7r/d, resulting in a primitive unit 
cell with the same volume as the original 3D Brillouin 
zone. We call the resulting unit cell the valley scattering 
unit cell. We define coordinates for each valley r so that 
it lies within the VSC. 

QW structures induce characteristically anisotropic 
scattering potentials, defined in the a-basis as V s (a, &, c), 
which can result in anisotropic intervalley scattering. 
Scatterers such as a single monolayer step in the side- 
wall of a QW or a miscut substrate with periodic side- 
wall steps are sharp on the order of a single lattice period 
and could in principle induce large momentum intervalley 
scattering in the A-space plane parallel to the QW. 

Recalling that the intervalley scattering vector is Q rr , 
we calculate the intervalley scattering matrix element, 

v T - = (Vv(r)|V s (a,6,c)|Vy(r)) (44) 
= <<Kr)tv(r)|F s (a, b, c)e l ^* +b+c ^(r) Uci r(r)) 

d 3 r^{r)u*{r)V s (a, b, c )e lQ "- c e lCl ^^ a+h) 



(r)u q -(r) 



(45) 



To emphasize the importance of the out-of-plane inter- 
valley scattering vector, we write Q rr in the exponential 
as the sum of an out-of-plane component Q^ r and the 
in-plane vector Qj[ r - If we use the standard planar 2D 
Brillouin zone view, we lose the representational impor- 
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FIG. 2: Depiction of the valley degeneracies for various QW orientations. Each row is labeled with its Miller index. Left 
column: 3D representation of the X-valleys showing both the CCC k x ,k y , k z and the transport basis k a , k c . In all diagrams, 
the doubly-degenerate X x ,y valleys are red, and the singly-degenerate X z valley is blue. For the (lll)-oriented AlAs QW, the 
three degenerate X x ,y,z valleys are purple. We note that the red and purple valleys exhibit robust valley degeneracy condition 
whilst the intersection point of the red and blue traces identifies the crossover degeneracy condition. The purple valleys satisfy 
the robust valley degeneracy condition as well as the crossover degeneracy condition trivially for all well widths. Center column: 
2D projection and 3D representation of the valley scattering unit cell described in the text. Right column: Valley degenerate 
ground energy as a function of QW width for a QW barrier of Alo.45Gao.55As at T = 4 K. Dashed lines represent the calculated 
strain energy shift of the respective electron valley relative to unstrained AlAs. The results of the heterostructure simulation 
software nextnano are depicted with scatter plots, and are verified with the analytical calculations shown as continuous lines. 
The (111) and (411) calculations assume that the piezoelectric field has been canceled by a top gate to result in a flat QW. 



tance of the out-of-plane scattering component. There- 



fore, we choose to depict the VSC in all three dimensions. 
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Specific examples will be shown in the following section 
which demonstrate how all possible intervalley scattering 
events can be easily identified in such a cell. 

VII. VALLEY DEGENERACY AND THE 
VALLEY SCATTERING UNIT CELL FOR ALAS 
QWS 

We apply the above analysis to the AlAs multi-valley 
system, determining the valley vectors, mass tensors, and 
strain tensors for the various growth directions, and then 
we calculate the valley subband energies as a function 
of well width W for each orientation. The three de- 
generate conduction band electron valleys are composed 
of six half valleys located at the A-points of the Bril- 
louin zone edge, with valley vectors q 100 = (2% /a) x, 
q 010 = (2tt/o) y, and q 001 = (2n/a) z for the X x , X y , and 
X z points, respectively. The superscripts express r as 
momentum-space directions after the notation of Van de 
Wallet. The AlAs electron valleys have anisotropic elec- 
tron mass, with heavy longitudinal mass mi = 1.1 m e and 
light transverse mass m t = 0.20 m c , respectively ^ 12 ' 13 . 
In the conventional crystal x-basis, the mass tensor m r 
is diagonal with = mi for the mass component par- 
allel to q T and mj^ = m t for the two transverse mass 
components repectively. 

Quantum confinement is created by sandwiching 
the AlAs QWs between aluminum gallium arsenide 
(Al^Gai-^As) layers which have a high aluminium con- 
tent x > 0.4. Properties of the barrier alloy are deter- 
mined by interpolating between AlAs and GaAs, either 
linearly or with a bowing term where applicable. We 
assume the reciprocal mass tensor in the Al^Gai-^As 
barrier layer to follow a linear interpolation for use in 

Eq. (dDj and jnjfsa, 

w Al x Ga (1 _ x) As = 0) w AlAs + i 1 - 2; ) w GaAs ( 46 ) 

For Al^Gai-^As, E^ is given by the relation 52 

^AUGax-.As =xEl iAs + {l-x)E£ aAs -bx(l-x) . (47) 

For Al K Gai_ x As 52 , the bowing term 5 ^ b T is 0.055 eV for 
i?Q aAs defined to be zero and E AlAs = 0.259 eV where r 
can be X x , X y or X z respectively. 

For all strain calculations, we assume AlAs to 
be strained with respect to the GaAs substrate because in 
typical structures the thickness of the intervening buffer 
Al^Gai^As layer is below the critical strain relaxation 
thickness of 0.5 /i m 58 ' 59 . Using Eq. (|2"2"|) we can calculate 
the in-plane strain from the lattice mismatch between the 
AlAs QW and the GaAs substrate 

, OGaAs — OAlAs /.„•> 

e,, - , (48) 

" OAlAs 

where ac a As = 0.564177 nm (0.565325 nm) is the lattice 
constant of the GaAs substrate and oaias = 0.565252 nm 



(0.566110 nm) is the lattice constant of the AlAs at K 
(300 K) 52 . The perpendicular strain component e' ± is 
calculated using Eq. (f23|) where Dq is the biaxial Pois- 
son ratio in the respective growth direction. The strain 
energy is determined from the absolute deformation po- 
tential for uniaxial strain at the X-point, = 6.11 eV 
in AlAs. 

In the following subsections, we compare and contrast 
the calculated results for the (001)-, (110)-, (111)-, and 
(411)-oriented QWs, as summarized in Fig. [5] In the left 
panel of the figure, the transport basis vector c is taken 
along the growth direction. The differently colored el- 
lipsoids distinguish the triple (purple), double (red), and 
single (blue) degenerate valleys, each of these satisfy the 
robust degeneracy condition given by Eqs. (|37|) or (|38l) . 
For the (001), (110) and (411) orientations, the strain en- 
ergy breaks the three fold valley degeneracy into twofold 
valleys and a single valley. The growth direction governs 
which of these valleys will have the lowest energy for 
a given range of well widths. For the (111) orientation, 
symmetry dictates threefold degenerate X x>y>z valleys for 
all well widths. 

The right panel of Fig. [2] shows the results of the an- 
alytical simulation of Eq which gives the valley subband 
energies as a function of the well width. The horizon- 
tal dashed lines indicate the pure strain splitting energy 
AE T , equivalent to the energy splitting in the asymptotic 
wide- well limit. The analytical results are confirmed by 
the semiconductor heterostructure simulation software 
nextnanc 60,61 . For the crossover QW widths Wq in- 
dicated for the (001), (110) and (411) orientations, all 
valleys satisfy the degeneracy condition of Eq. (|3"3")l simul- 
taneously. The intervalley scattering momentum vectors 
Q TT for all the facets are listed in Table HI1 



A. (001) AlAs 

For the (OOl)-oriented QW, the crystal axes are iden- 
tical to the growth axes, Fig. [2] (top row), thus the co- 
ordinate transformation matrix is the identity matrix 
R 001 — I. The inverse mass tensor follows trivially. The 
various mass parameters listed in Table Hi] are then used 
to calculate the respective energy terms in Eq. (0. The 
biaxial Poisson ratio derived from Section IIIII is given 
by22 Dff = 2^jj, and the diagonal strain matrix has 
D® 01 = D® 01 = 0. The strain tensor in the a-basis is 

"10 



-0.854 

The results of the Section HTT1 analysis for (OOl)-oriented 
QWs including 2D density-of-states and cyclotron mass 
for each valley are provided in Table |TTJ 
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TABLE II: Reciprocal mass tensors w' T for the various valleys in the transport basis for (001)-, (110)-, (111)-, and (411)- 
oriented AlAs QWs. Intervalley scattering vectors Q TT are listed, as are the 2D density-of-states (cyclotron) masses relative to 
the free electron mass. The 2D density-of-states values n^D are combined for robust-degenerate valleys. 



The center panel of Fig.[2](001) depicts the VSC and its 
2D projection. The 2D sublattice of coplanar reciprocal 
lattice points is square symmetric, with the X z valleys 
lying at the corners of the unit cell and the X XtV valleys at 
the edges. The resulting VSC illustrates that all valleys 
can be connected by purely in-plane scattering events, 
which is not obvious from the usual depiction of (001) 
valleys. This can be verified in Table |TT] where the c- 
component of the scattering vector QY is zero for all 
intervalley scattering vectors. 



As seen in Fig. [2] (001, right), strain alone as indicated 
by the horizontal dashed lines contributes an energy shift 
of AE 100 > 010 = -11.62 meV for the doubly-degenerate 
X XtV valleys while raising the singly-degenerate strained 
X z ' valley by AE° 01 = +9.93 meV, a net difference of 
A = 21.55 meV. The ordering of the degeneracy changes 
at the crossover well width Wq 01 = 5.4 nm. 



Because (001) biaxial strain induces no shear compo- 
nent in the crystal basis, there are no piezoelectric fields 
for this facet. 



B. (110) AlAs 

For the (HO)-oriented QW in Fig. CD (110, left), we ap- 
ply the rules of axis identification described in Section 
IIII1 a lies along the lowest Miller index [001] direction, 
making the X z valley in-plane, and b along [110] com- 
pletes the right handed coordinate system. The coordi- 
nate transformation matrix R 110 is 



R 



no 



1 
-4= =1 

J- ^- 



(50) 



The biaxial Poisson ratio derived from Section|TTT]is given 
by^ 2 . D¥ = Cll + 3ci2 ~ 2c44 , and the diagonal strain matrix 

J U Cii+Ci2+2c 4 4 ' o 

has D\ w = D\ w = 0. 

The center panel of Fig.[2] (110) depicts the VSC and 
its 2D projection. The strain tensor is diagonal in the 
a-basis and given by 



1 
1 








-0.616 



(51) 
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The results of Section HTT1 analysis for (llO)-oriented QWs 
including 2D density-of-states and cyclotron mass for 
each valley are provided in Table HH 

The 2D sublattice which defines the VSC has a 
centered-rectangular symmetry. The X z valley lies in 
this plane and the X XiV valleys lie outside of the central 
plane. In the VSC, it is clear that whereas the X xy val- 
leys can scatter amongst each other with purely in-plane 
scattering, the X z valley is isolated and requires an out- 
of-plane component for intcrvalley scattering. Table |TT] 



a-basis is 



reflects this result since only the Q 



100,010 



scattering vec- 



tor has a zero c-componcnt. 

As seen in Fig.[2](110, right), the singly-degenerate X z 
valley has the lowest strain energy E 001 — — 12.94 meV 
whereas the doubly-degenerate strained X XtV valleys have 
an energy of E 100 ' 010 = — 3. 55 me V due to the smaller 
compressive strain. A smaller strain differential of A = 
10.77 meV is observed between the valleys in this orien- 
tation as compared to the (001) orientation. We observe 
the valley degeneracy crossover at Wg 110 = 5.3 nm. 

A purely in-plane piezoelectric field arises in (110) 
QWs due to the nonzero shear component of the strain 
tensor in the x-basis. When Eq. (|23|) is transformed to 
the x-basis, the shear components take the fori»S& 



p — Cll— 2C12 J 

tx V Cn+Ci2+2C 4 4 t || 



(52) 



leading to an electric field in the QW plane parallel to 
the [100] a-direction 



E' n = 



£s^0 



(53) 



This electric field £' a = 7.942 x 10 6 V/m plays no role 
in quantum confinement so that the flat QW assumption 
remains valid, and the field is screened in-plane by the 
electrons in the QW within a Thomas-Fermi screening 
length of the sample edges at ±a. 



C. (Ill) AlAs 

For the (lll)-oriented QW, a is chosen along the lowest 
Miller index [110] direction perpendicular to the growth 
axis c, and none of the valleys lie within the plane, b 
is chosen along the [112] direction to complete the right 
handed coordinate system^ .. The coordinate transfor- 
mation matrix R 111 is 



R 



J_ =1 n 

11-2 

Ve Ve Ve 
ill 

\/3 V3 V3 



(54) 



The biaxial Poisson ratio derived from Section|TTT]is given 
by^ 2 . D¥ = 2 £n±|a2^|a4 d th t in matrix is di- 

J U Cii+2ci2+4c 4 4 ' 

agonal with D\ u = D^ 11 — 0. The strain tensor in the 



1 
1 








-0.55 



(55) 



The results of the Section ITO1 analysis for (lll)-oriented 
QWs including 2D density-of-states and cyclotron mass 
for each valley are provided in Table IPI1 The (111) QW 
valleys remain threefold degenerate due to equal contri- 
butions of the strain and confinement energies for all elec- 
tron valleys as seen in Fig. [U (111) left. 

The center panel of Fig.[2](lll) depicts the VSC and its 
2D projection. The planar sublattice is hexagonal, with 
six half-valleys located at the center of the hexagonal 
facets. The resulting VSC shows that the valleys are all 
connected with coplanar scattering vectors, as seen in Ta- 
ble [IT]from the zero c-component of all Q rT . Fig. [2] (111, 
right), shows that all the valleys in this triply-degenerate 
system are equally strained to AE 1100 ' 010,001 = — 7.0meV 
and the robust degeneracy condition given by Eq. (|33)) is 
trivially satisfied for all well widths. 

A purely out-of-plane piezoelectric field arises in (111) 
QWs due to the nonzero shear component of the strain 
tensor in the x-basis. When Eq. (|2U)) is transformed to 
the x-basis, the shear components take the for m 24 i 27 



-Gii — 2Ci2 



C u + 2C 12 + 4C 44 



"II 



(56) 



leading to an electric field perpendicular to the QW plane 
parallel to the c-growth axis: 



£' = 



-2e xA \/ 3 



(57) 



We will assume that this electric field £' c — 8.794 x 
10 6 V/m will be canceled by a gate voltage, so that the 
flat QW assumption remains valid. 



D. (411) AlAs 

To demonstrate the utility of the formalism described 
in this paper, we now extend our analysis to the less con- 
ventional (411) AlAs QW which has nonetheless shown 
important results in the liter ature^S 2 '. Applying the same 
rules of axis identification, a is chosen along the [OlT] di- 
rection which is the lowest Miller index along the plane 
of the QW. b is chosen along the [122] direction to com- 
plete the right-handed coordinate system. The coordi- 
nate tranformation matrix R 411 is 



R 



411 



1 

V2 
2 
3 
1 

SB 



o =k 



3 
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(58) 
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FIG. 3: Illustration of shear strain for the (411)-oriented QWs 
due to the nonzero strain-to-shear ratio Df 11 . Ao, Ab and Ac 
represent strain displacements along the respective vectors of 
the transport a-basis, and a represents the shear vector in 
the plane of the QW, which for the (411)-oriented QWs is 
parallel to the b-axis. 



Using the derivation for the strain ratios from Section lllTl 
for the (411)-oriented QWs, we obtain Dq 11 = 0.775 and 
D^ 11 = 0.176. The explicit relations for the biaxial Pois- 
son ratios for AlAs (411)-QWs are 



D, 



411 



6(cn + 2ci 2 )(4cii - 4ci 2 + 19c 44 ) 

8cfj - (16ci2 - C44)(ci2 + 2c 44 ) + cn(8ci2 + 145c 44 ) 



D 



411 



-15\/2(eii + 2ci 2 )(cn - C12 - 2c 44 

8^ 



(16ci2 - c 44 )(ci2 + 2c 44 ) + cn(8ci2 + 145c 44 ) 



(59) 



The results of Section ITU analysis for (411) including 
2D density-of-states and cyclotron mass for each valley 
are provided in Tableim For this high-index facet, the 
strain tensor is nondiagonal in the a-basis and given by 



1 







1 0.176 
0.176 - 0.775 



(60) 



The biaxial strain in the structure grown along [411] 
gives rise to nondiagonal components in the strain tensor. 
The effect of shear is depicted in Fig. [3] by the shear 
vector a (Eq. (|32|l) which lies in the plane of the QW and 
is parallel to the 6-direction. The degeneracy condition 
takes into account the shear component Df 11 and is given 
by Eq. (01D. 

Fig. H] (411, center) depicts the VSC and its 2D pro- 
jection. The 2D planar sublattice in Fig. [5] (411, center) 
has centered rectangular symmetry. 



From Fig. [2] (411, right), we obtain a crossover width 
Wq 11 = 5.0 nm, and the valley occupation as a func- 
tion of well width follows an analogous discussion as for 
the (110) case. However, we observe that the X x val- 
ley is singly-degenerate and the X V:Z valleys are doubly- 
degenerate. The strain component of the X x valley en- 
ergy is AE 100 = 6.82 meV, and that of the X y . z val- 
leys which satisfy the robust valley degeneracy condi- 
tion is AE 010 = -12.29 meV. A strain differential of 
A = 19.11 meV is obtained. The two degenerate val- 
leys have coplanar momentum scattering as seen in the 
vector Q 010 ' 001 which has zero c-component. 

A piezoelectric field with both in-plane and out-of- 
plane components arises in (411) QWs due to the nonzero 
shear component of the strain tensor in the x-basis. 
When Eq. (f2"3"| is transformed to the x-basis, the shear 
components take the formal 



W 1 



li 



-(Cii+2Ci 2 )(17Cii-17Ci2+2C 4 4) 

(16Ci 2 -C44)(C 12 +2C44)+C 11 (8Ci 2 + 145C 4 4) II 



(61) 



-(Cii+2Ci 2 )(8Ci 



-8Ci 



-7C 4 4) 



8Cjf 1 -(16Cl2-C44)(Cl2+2C44)+Cll(8Cl2+145C44) II 



leading to an electric field with components both in the 
QW a-b plane as well as parallel to the c-growth axis: 



£' = 



2e 



■r, 4 



£s£o 



~ 3 c xy 



V2 €zX 

3^ 



(62) 



We will assume that the out-of-plane piezoelectric field 
component £' c = 1.783 x 10 6 V/m will be cancelled by a 
gate voltage so that the square QW assumption remains 
valid. The piezoelectric field component £' a is zero, and 
the in-plane component £' b = 3.761 x 10 6 V/m will be 
screened at the sample edges by the electrons in the QW. 



VIII. THE VALLEY SCATTERING UNIT CELL 
AND ANISOTROPIC INTERVALLEY 
SCATTERING FOR ALAS QWS 

This section first discusses the differences between the 
VSC and standard depiction of the 2D Brillouin zone, 
and then determines the inter- to intravalley scattering 
ratio for valleys near crossover degeneracy. In the VSC 
description for the case of (lll)AlAs, Fig. [4] shows the 
(lll)AlAs VSC, as well as the 2D Bravais lattice of these 
cells. For comparison with this hexagonal unit cell, the 
standard 2D Brillouin zone is plotted as a smaller grey 
hexagon within this latticed. The hexagons A, B and C 
are intended to illustrate the relative positions of three 
identical planar hexagonal lattices, displaced from each 
other. When these three layers are stacked in sequence 
ABCABC out of the plane we recover the full 3D recip- 
rocal lattice. Whereas the standard 2D Brillouin zone is 
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(110) k 




FIG. 4: 3D representation of the valley scattering cell for (lll)AlAs as well as 2D lattice of such cells. In AlAs, the valleys are 
located at the q = (2n/a)( x, y, z) The standard depiction of 2D Brillouin zone (small shaded hexagon)-i projects together 
three identical but laterally and vertically displaced planar hexagonal lattices, translated according to the representative 
hexagonal cells A, B and C. The area of standard 2D Brillouin zone is one third that of the valley scattering cell since its height 
is three times that of the valley scattering cell. Arrows illustrate the planar intervalley scattering vectors Q Tr for (lll)AlAs. 



defined from the overlay of all three layers, we define the 
VSC from the 2D Brillouin zone of a single layer only. 
The area of the hexagonal VSC for (111) is thus three 
times larger than that of the traditional hexagonal unit 
cell. The coplanar intervalley momentum scattering vec- 
tors are shown with arrows in the VSC. 

The VSC is helpful in depicting intervalley scatter- 
ing events, so we will discuss the two main sources of 
single-particle scattering for AlAs QWs, namely inter- 
face roughness and alloy disorder scattering in the bar- 
rier. We find below that near valley degeneracy, interface 
roughness scattering does not result in any significant in- 
tervalley scattering, and in the presence of alloy scatter- 
ing in the barrier walls, the scattering rate for crossover 
intervalley scattering is significantly repressed relative to 
robust intervalley scattering, and this suppression factor 
is calculated. 

We assume low enough temperatures such that acous- 
tic and optical phonons can be neglected. Elastic in- 
tervalley scattering can only occur if the valley splitting 
energy is less than the Fermi energy at low temperatures, 
\E T — E T \ < Ep, or the thermal energy at high tempera- 
tures \E T — E T \ < k^T. Two independent processes can 
in principle result in elastic intervalley scattering: inter- 
face roughness scattering and alloy disorder scattering 
due to wavefunction penetration into the barrier alloy. 

Interface roughness leads to local changes in the well 
width A(r) at an in-plane position r = (a, b), and the 
resulting fluctuation in the quantization energy causes 
a scattering potential^. The roughness is characterized 
by the rms average displacement of the interface A and 



the roughness correlation length along the interface plane 
A, and is expressed as the Gaussian autocorrelation 
function^ 



< A(r)A(r') >= A 2 exp - 



\r-r'Y 



(63) 



where < ... > means an ensemble average. The Fourier 
transform of the roughness autocorrelation at the inter- 
valley scattering vector Q rr is 



<| Aqt- | 2 >=^(A 2 A 2 )exp 



(AQ TT ) 2 



(64) 



Following Quang et al£^& the scattering potential U-rs 
for the QW with a square potential and a well width W 
is given by 



< 



'RS 



: >=A 



1,1 W,cc 



2W 3 



<l A 



(65) 



where w' ^ cc is the inverse confinement mass in the well 
perpendicular to the interface and A is a unitless nor- 
malization constant. Under optimal epitaxial conditions 
of slip-step growth, the AlAs roughness disorder poten- 
tial will have correlation lengths of tens of lattice periods 
A = Na. Since the intervalley scattering vector Q rr is of 
order a Brillouin zone boundary it /a and the intravalley 
scattering vector is of order zero, the ratio of intervalley 
scattering to intravalley scattering is of order exp(— V 2 ), 
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which even for a conservative estimate of N = 5 lattice 
periods of autocorrelation already yields a suppression of 
more than 10 orders of magnitude for intervalley scat- 
tering. Thus intervalley surface-roughness scattering is 
significantly suppressed with respect to intravalley scat- 
tering. 

Alloy scattering is the other possible mechanism for 
elastic intervalley scattering. The standard model for 
alloy disorder scatterin g 65 ! 67 with variable alloy compo- 
sition x{z) is described by the autocorrelation function 



<l un \ 2 > 



n 3 V 2 
"0 v 



lp 2 T (c)x(c)[l-x(c)]Tj4(c)dc 



(66) 

where clq is the in-plane lattice constant, and Vo is the 
spatial average of the fluctuating alloy potential over the 
alloy unit cell. In AlAs QWs, only tails of the wavefunc- 
tion have an alloy content in the Al x Gai_ x As barrier. 
One notes that the wave functions of electrons in valleys 
with a light confinement mass penetrate to a greater ex- 
tent in the barrier and will have a significantly larger alloy 
scattering potential as compared to the heavy confine- 
ment mass. Assuming valley degeneracy for all valleys 
and plane wave solution for the electron wavefunction in 
the QW barrier, the scattering matrix element can be 
written as 



<\U 



AD 
n 3 V 2 



> = 



A 2 -2^c 



x(c)[l - x(c)]A±e 



2 „-2k t c 



dc 



,(67) 



where A T = \^ T (W/2)\ and A- = \iPr(W/2)\ arc the 
boundary values of the wavefunction at the barrier, ob- 
tained by solving for boundary conditions of the square 
wave potential. 

We get a quantitative estimate of the ratio of the al- 
loy scattering 7 X between crossover-degenerate valleys 
labelled t\ and T3 compared to the scattering 711 be- 
tween robust-degenerate valleys labeled t\ and T2 in the 
Al a Gai_ x As barrier by calculating 



<l U 



AD 



> 



<l Uj^ I 2 



_ = 1* = 

> 7|| 

fw A 2 T1 e" 2 ^ 1 c x(c) [1 - x(c)] A 2 T3 er 2 ^ c dc 
^ lc x(c)[l - x{c)]A 2 e- 2KT2c dc 



(68) 



f A 2 

Jig- A TX 



where c is the direction of integration out-of-plane. For 
the case of (001)-, (110)-, and (411)-oriented AlAs QWs, 
2* = 0.0012, 0.0691 and 0.0235 respectively. The val- 
ley with the larger confinement mass penetrates less into 
the barrier region and results in a smaller scattering in- 
tegral ratio. The VSC makes it easy to visually identify 
this anisotropic scattering possibility along the various 
orientations. 

Electron-electron intervalley scattering is outside the 
scope of this article and will be treated elsewhere. It 



is expected to be small since the Fermi wavevector kp 
is much smaller than the intervalley scattering vector— 
Q rT . 



IX. VALLEY DEGENERACY IN MISCUT 
SAMPLES 

In this final section, we show how the analysis of Sec- 
tion IIIII can be applied to miscut samples to determine 
the projected masses, ground state energies, and strain 
energies. Miscut samples are prevalent in the literature, 
for example (111) Ga As/ AlAs growth has been shown to 
have superior morphology with a miscut angle from 0.5° 
to 4° and in general, intentional miscuts improve growth 
quality due to slip-step growth 69-72 . In Si, recently in- 
vestigated hydrogen terminated (lll)Si miscut surfaces 
have shown very high mobility^ and the wafer miscut is 
expected to break the valley degenerac y 43 ! 74 . 

Miscut samples are characterized by two angles. One 
angle (f> designates the azimuthal angle in the a-b plane 
relative to the a-axis towards which the plane is tilted, of- 
ten expressed as an in-plane Miller index tilt direction T. 
The other angle 9 designates the polar tilt angle relative 
to the surface normal c. Experimentally, the angles can 
be deduced from atomic force microscope images of sur- 
face monolayer steps, where <j> is oriented perpendicular 
to the steps, and 9 = atan(a/2L) is deduced from know- 
ing the monolayer thickness a/2 and the average width 
of the monolayer steps L. We define a new coordinate 
transformation matrix R M * for the miscut samples, by 
introducing the azimuthal and polar rotations R^ and 
Re, respectively: 



R-T — Rrf 



R» 



cos(/> — sin</> 
sin0 cos(p 
1 



cos# sin# 

1 
— sin(9 cos6> 



(69) 



(70) 



We obtain a new coordinate transformation matrix 
for calculating the transport parameters and the 
strain tensor analysis of miscut samples 



R-m — (R^RflR^R 



M 



(71) 



Note that a VSC cannot be defined for miscut samples, 
but only for those with a pure Miller index. 

We outline here the procedure for determining valley 
subband energies under miscut. First, the transforma- 
tion matrix Rg ^ from Eq. (|71[) determines the inverse 
mass tensor of each valley with Eq. ©, and the out- 
of-plane inverse mass tensor component can be used in 
Eq. (|10[) to determine the ground confinement energies 
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of each valley. The in-plane inverse mass tensor is then 
calculated with Eq. (fT4")) and used in Eq. (fT5l) to obtain 
the kinetic energy. The stress tensor can be determined 
by applying ± in Eq. (T25|) . and the strain ratios Df 1 
can be determined by minimizing the elastic energy with 
Eqs. ([28j) - ([30l) . Finally, the strain energy shift can be 
deduced from Eq. (|19|) and added to the ground quan- 
tum confinement energies and kinetic energies to arrive 
at the final valley energies. Any piezoelectric fields can 
be determined by transforming the strain tensor to the 
unprimed frame e and inserting the appropriate shear 
components into Eq. (|4"2"j) . 



X. CONCLUSION 

In conclusion, we review the advantages and limita- 
tions of the formalism we have developed for valley sub- 
band energy calculations. We first derived how the valley 
index is a valid pseudospin index in a multi-valley sys- 
tem. A key element of our model is the definition of two 
relevant bases: the conventional crystal x- and the trans- 
port a-basis. We use this notation to find the projected 
in-plane and out-of-plane effective masses, the electron 
subband energy and degeneracy, piezoelectric fields, and 
the scattering vectors. 

There are five competing energy scales in multi-valley 
QWs namely in-plane kinetic energy, confinement energy, 
strain energy, piezoelectric energy, and Fermi energy, all 
of which influence valley occupancy. Strain in the QW 
breaks valley degeneracy and results in an energy differ- 
ential, AE T of - 20 meV between non-degenerate val- 
leys. The strain energy is independent of the QW width 
unlike the quantum confinement energy which increases 
as the QW width decreases and is inversely proportional 
to the confinement mass; typical confinement energies Eq 
are of the order of 20 meV for a narrow QW about 5 nm 
wide and 5 meV for a 20 nm wide QW. We defined a ro- 
bust valley degeneracy condition to identify valleys which 
are degenerate independent of well width, and a crossover 
degeneracy condition where valleys are degenerate only 
for a particular crossover well width Wo- We showed that 
piezoelectric fields in the QW and barrier materials can 
arise in certain facets such as (111) and (411), and cause 
piezoelectric energy shifts E^ z as large as 50 meV for a 
(111) QW about 10 nm wide, and that such a piezoelec- 
tric field can be canceled by an external gate voltage. 
The smallest energy scale in these QWs is the Fermi en- 
ergy which is a few meV depending on the degeneracy 
for electron densities around U2D ~ 3 x 10 11 cm~ 2 and 
assuming an effective mass mjrj ~ 0.3m e . With such a 
small Fermi energy, valley occupation is strongly defined 
by the interplay of the various other energy scales in the 
problem. We once again emphasize that it is due to the 
similarity of all the energy scales that none of them can 
be ignored and must be treated carefully to calculate val- 
ley subband energy. 

Valley degenerate systems have an extra scattering 



channel not present in single valley systems, namely 
inter-valley scattering. We identified the main inter- 
valley scattering mechanisms as the interface roughness 
and alloy disorder scattering due to wavefunction pene- 
tration into the barrier and calculated the suppression 
factor at the crossover well width Wo for inter-valley 
scattering between different sets of valleys. We also 
defined a valley scattering primitive unit cell to easily 
identify scattering events between the robust-degenerate 
and crossover-degenerate electron valleys. We drew the 
VSC for the special case of AlAs grown along four differ- 
ent orientations, the high-symmetry facets (001), (110), 
and (111), as well as a low-symmetry facet (411) to 
demonstrate the utility of our model. Furthermore, we 
explained the relevance of the VSC description for vi- 
sual identification of anisotropic inter-valley scattering 
in AlAs QWs. 



In the final section we demonstrated the power of our 
formalism to determine valley degeneracy for arbitrary 
substrate orientation, without being restricted to Miller 
indexed planes as is prevalent in the literature. We con- 
sider the example of a miscut sample and define the addi- 
tional notation required to address angle deviation from 
conventional growth axes. We then detail the procedure 
to calculate all the ground energy parameters as well as 
the strain tensor analysis for miscut samples. It is worth 
repeating that the formalism developed here is sufficient 
to calculate valley degeneracy for any substrate orien- 
tation, and need not be aligned with any Miller index 
facet. 



There are some important limitations of our formal- 
ism. Firstly, we assume layer thicknesses are all within 
the strain relaxation limit such that there are no strain 
induced defects below or within the QW. These would 
relax the lattice constant relative to the substrate, thus 
one would have to first determine the adjusted lattice 
constant at the QW layer and then apply the formalism 
developed here to deduce valley degeneracy. Secondly, 
we assume nearly empty QWs and do not calculate the 
Schrodinger equation self-consistently with the Poisson 
equation. However our formalism for calculating the var- 
ious valley subband energies can be applied to standard 
self-consistent solvers to obtain the proper Hartree solu- 
tion. Lastly, whereas the formalism takes into account 
the change in the confinement potential due to piezoelec- 
tric fields, it does not address larger structural consider- 
ations such as modulation doping layers, surface pinning 
potentials, electrostatic charges, and piezoelectric fields 
in the barrier, all of which would have to be calculated 
together to determine the confinement well potential and 
Fermi energy. We note, however, that the piezoelectric 
equations provided here can be used to determine the 
piezoelectric fields in the strained barriers. 
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Appendix: Analytical equations of the strain tensor 
and the strain ratios for arbitrary substrate 
orientations 

In Section [III1 we derive the strain ratios Df 1 by min- 
imizing the free energy of the strained layer. This sec- 



tion presents a generalized treatment to express these 
strain ratios for an arbitrary substrate orientation with 
cubic symmetry, given its coordinate transformation ma- 
trix R , the elastic constants Cijki in the CCC x-basis 
and the in-plane strain ej| in the transport a-basis. 

We first obtain the fourth order rotated elastic stiff- 
ness tensor c' mnop in the a-basis by rotating the fourth- 
rank elastic stiffness tensor c^ki from the CCC basis 



x 



Cijkl 

to the transport basis a 



(a, 6, c) 



ijkl 



M i 
mi 



(A.l) 



where R is the coordinate transformation matrix de- 
fined in Eq. (Q. The rotated fourth-rank tensor c' mnop 
in the CCC x-basis should then be mapped to using 
the Voigt notation. The strain tensor components in 
the transport a-basis are given by 
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where the denominators in Eq. (IA.4I) . Eq. (|A.5|) and 
Eq. (|A.7[) are always nonzero, the coefficients a, /? and 7 
are first order, and the coefficients A, k, rj, lo and fj, are 
second order in the elastic stiffness tensor components 



a 


— — (^13 + ^23) 


(A.7) 




— ~ (^14 + C24) 


(A.8) 


7 


= -ic' 15 + c 25 ) 


(A.9) 


A 


— 2C 33 C 44 


(A.10) 


K 


= 2C" 33 C' 55 


(A.11) 


V 


= 2 (C 33 C 45 — C 34 C 35 ) 


(A.12) 


LO 


= - C' 3i a 


(A.13) 


/' 


= C' 33 -i-C' 3b a . 


(A.14) 



